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A Lagrangian analysis explains the numerical results of Valentini and Westman 
(2005) which demonstrate that an initially arbitrary particle density, stirred by the 
field of de Broglie velocities associated with a Schrodinger wave function, relaxes to 
the Born probability density provided both densities are coarse-grained. 
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1. Introduction 



The de Broglie-Bohm interpretation of quantum mechanics (Bell, 1987; Bohm and 
Hiley, 1993; Holland, 1993) continues to attract critical appraisal (e.g., Passon, 2006) 
in spite of being beset by several difficulties. The most immediately apparent is that 
the initial density for ontological particles may not be freely chosen, but must be the 
Born probability density at that time if the two densities are to agree subsequently. 
The resolution of the difficulty is indicated by numerical calculations of Valentini 
and Westman (2005) who, using a Schrodinger wave function demonstrate that an 
arbitrary initial choice for the ontological or Bohm particle density p relaxes to the Born 
probability density provided both densities are coarse-grained. The relaxation owes 
to the stirring action of the velocity field v = v[x, t] prescribed by de Broglie in terms 
of ip: 

v=3/\i^\\ (1) 

where j is the Schrodinger current: 

3 = (n/m)Im(V'*VV') . (2) 

The demonstration by Valentini and Westman (2005) of the relaxation of the coarse- 
grained particle density to 'quantum thermal equilibrium' is a major development for the 
de Broglie-Bohm ('dBB') theory. Valentini (2008) has proposed that particles emitted 
in the early universe may not yet be in equilibrium. The analogy of the dBB formalism 
with fluid dynamics suggest borrowing the Lagrangian representation from the latter 
(e.g., Holland, 2009). It is shown here in §2, using the Lagrangian representation, that 
the relaxation demonstrated by the numerical computations of Valentini and Westman 
(2005) is in general to be expected. 



2. Relaxation: Schrodinger equation 

2.1. particle kinematics 

The Eulerian field v — v[x,t\ of the de Broglie velocity (1) defines particle paths 
X — P{a, s; t) through the Lagrangian prescription 
dP 

— {a,s;t)=v[P{a,s;t),t], (3) 
subject to the general initial condition 

P(a, s; s) = a . (4) 

Note that both t > s and t < s are allowed. In particular, if a; = P{a,s;t) then 

a — P{x,t]s) . In the interest of clarity, Eulerian arguments arc denoted by [x.t], 
Lagrangian arguments by (a, s; t) . The Jacobian of the transformation a P{o,, s; t) 
is denoted by J = J(a, s; t) where 

d{a\a^,a^) ' ^ ^ 
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In particular, J(a, s]s) — 1 . It is a kinematic identity (e.g., Lamb, 1932; Bennett, 
2006) that 

^{a, s; t) = J{a, s; t)V ■ v[P{a, s; t),t] . (6) 

Arguments will not always be displayed with such care subsequently, but the meaning of 
the partial derivatives may be inferred from context. For any field t] , the following 
relationship expresses the rate of change of F following the motion both in Lagrangian 
and Eulerian terms: 

OF dF 

— {a, s; t) = -Q^[P{ci, s; t),t] + v[P{a, s; t),t]- VF[P{a, s; t),t] . (7) 

2.2. Bohm particle density 

According to the dBB theory, the particle density p — p[x, t] obeys the Eulerian 
conservation law 

^ + V-M = 0. (8) 
This may be rearranged as 

^ + v-Vp=-pV -v , (9) 
which by virtue of (6) and (7) has the Lagrangian form (Lamb, 1932, Art. 14) 

^{a,s-t) = 0. (10) 

2.3. Born probability density 

All solutions ip = ip[x,t] of the Schrodinger equation 

where V is the potential field (e.g., Feynman et al , 1965), conform to the dynamical 
identity 

^ + V.i = 0, (12) 

where the Schrodinger current j is defined in (2). This identity may be rearranged to 
resemble an Eulerian 'probability conservation law' involving the de Broglie velocity (1): 

^ + V.{v\m = 0, (13) 
which has the Lagrangian form 

^(....)^0. (14, 

It is clear from (10) and (14) that p and It/'P are identical at all times t if and only if 
they are identical at one time s, say s = . It is tacitly assumed that the particles 
are inside the domain D for all time, which is assured if i/j*'Vijj is real on the boundary. 
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This condition is satisfied by a variety of natural clioices of boundary conditions for the 
Schrodingcr equation (11). It also suffices that -0 be periodic in space (Valentini and 
Westman, 2005). 

2.4- Coarse-graining 

The stirring of the particle density by relatively smooth fields of de Broglie velocity leads 
in general to very fine structure in the density field (Valentini and Westman, 2005). The 
structure may be suppressed by spatial smoothing or 'coarse-graining' in space: 

{p[x,t]}^w-^ I p[y,t]dy, (15) 

J cell 

where the region of integration is a small but finite rectangular spatial cell centered on 
X , the infinitesimal volume element dy is dy^dy'^dy'^ and w is the cell volume. Each 
point y in the cell at time t lies on a path from a point b — P{y, t; 0) somewhere inside 
D at time t — , that is 

{p[x,t]} = vj-^ j p{b,0;t)J{b,0;t)db , (16) 

where p(6, 0; t) = p[P{b, 0; t),t] . The region of integration <S = S{x, t) is the pre-image, 
at time t = 0, of the cell centered on x at time t . In general 5 is a highly convoluted 
but simply connected region. Equation (9) shows that there is a tendency for particles 
to concentrate in zones of divergence (V ■ v > 0) as time decreases, but the zones are 
in general so transient that the stirring spreads the cell throughout D , albeit into a 
convoluted shape. By virtue of the conservation law (10), the coarse-grained density is 

{p[x,t]} = w-' J^p[b,0]db, (17) 

since p(6, 0; 0) = p[P(6, 0; 0), 0] = 0] , and J(6, 0; 0) = 1 . Now at any time t the 
particle density p is normalized over D, and in particular at t — , thus 

p[b,0]^W-' + p'[b,0], (18) 

where W is the volume of D , and the integral of p' over D vanishes. Substituting for 
p in the right hand side of (17) yields 

{p[x, t]} = {Ww)-' J db+... . (19) 

The ellipsis in (19) is small, owing to the integral over the 'tentacles' of S samphng 
and summing the signed values of p' . The integral in (19) is the volume of <S, denoted 
W = yV{x, t), thus the the coarse-grained particle density is found to be 

{p[x, t]} = {WwyWix, t] + ...^ W-^{J{x, t;0)} + ... . (20) 

By an identical argument it follows also that 

{\ij\^[x,t]} = {Ww)-'W[x,t] + ... = W-^{J{x,t;0)} + ... . (21) 

In general, It/'P is far less affected than p by coarse-graining (Valentini and Westman, 
2005). In any case (20) and (21) explain the relaxation of the coarse-grained Bohm 
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particle density, irrespective of its initial values, to the coarse-grained Born probability 
density. They further explain why 1-01^ may be used as a probability density. 

3. Discussion 

The elementary manipulations leading to (20) and (21) are not found in, for example, 
Tolman (1938). The density considered there is that of particles in phase space, where 
the flow is solenoidal owing to the underlying classical dynamics being Hamiltonian. 
Indeed, if x in (20) represents a point in phase space, and if J = 1 , then p[x, t] relaxes 
to the uniform density assumed by Tolman. The analysis in §2 does not depend upon 
the dimensionality of x , so relaxation also results from stirring by the current field of 
a Schrodingcr wave function on a higher-dimensional configuration space (xi, 332) • • •) • 
There are cases of interest in which relaxation does not occur, the most obvious being 
vanishing de Broglie velocity: v = 0, as in a standing Schrodinger wave. A uniform but 
time dependent flow would preserve cell shape. Rigid rotation: v — Vtxx, which would 
also preserve the cell shape, is not of interest since it is rotational: V x v — 2f2 . On 
the other hand a steady flow can stir effectively unless it restricts the cell particles to a 
region in which p' , the deviation of the initial particle density from its spatial mean, is 
one-signed. 

The stirring efficiency of the de Broglie velocity field controls the relaxation. The 
velocity field (1), associated with the Schrodinger wave function = exp{iS / h) , 
may also be expressed as v = m'^VS . The dynamics of this 'fluid' are revealed by 
decomposing the Schrodinger equation (11) for the amplitude and phase of ip (e.g., 
Holland 1993). The dynamics of V-S" resemble those of an inviscid fluid. However the 
held (1) is irrotational, and thus cannot be turbulent in the usual sense. The held (1) has 
divergence S — V ■ v , and the 'fluid dynamics' of this divergence may be rearranged as 
d6/dt+6'^ + .... = .... , where the time derivative follows the motion. The Riccati operator 
admits the possibility of blow-up in a finite time: 1/S{t) = 1/5(0) + t . On the other 
hand ip is the solution of the linear Schrodinger equation (11), and so the wavenumber 
spectrum of ijj is determined by the initial wave function and by the potential. Valentini 
and Westman (1995) report highly intricate and rapidly separating paths derived from 
a relatively smooth initial wave function and a uniform constant potential. 

Coarse-graining is crucial for relaxation of the Bohm particle density to quantum 
thermal equilibrium, and the formalism of a density at a single point is misleading. No 
matter how small the cell size, it will eventually (to speak in terms of time increasing) 
disperse into a space-exploring shape. A single particle remains forever a single particle. 
That is, relaxation is an effect of relative dispersion, not absolute. The relaxation time 
scale is not readily estimated in general. If the velocity field that stirs the particles 
is rich in fine structure then particles disperse independently, each taking a random 
walk, leading to root mean square particle separations growing like 0(t^/^) . If the 
dominant stirring effect is on relatively large scales, then relative dispersion owes to the 
large-scale rate of strain and separations grow exponentially in time. Random walks 
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through simple shear flows are ubiquitous in the ocean and atmosphere, with particle 
separations growing like 0(t^/^) . Random walks through a large-scale divergence lead 
to exponentially growing separations. In conclusion, the relaxation time scale cannot 
be expected to have a simple parameter dependence. 
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